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combustion  in  high-speed  flow  is  an  effective  mixing  of  gaseous  fuel  and  oxidant  providing  full  combustion  during  the  resident  time  in 
combustor  chamber.  Generally  speaking  in-flow  mixing  consists  of  two  main  mechanisms:  kinematics  mixing  (convection)  and  molecular 
diffusion.  The  latter  is  a  final  stage  of  mixing  needed  to  provide  extended  combustible  mixture.  Kinematics  mixing  is  more  effective  in  vortex 
flows;  the  higher  vorticity  the  better  kinematics  mixing  which  means  the  larger  fuel/oxidant  interface  surface.  The  second  mixing  mechanisms 
-  diffusion  forms  a  reacting  volume  estimated  as  the  interface  surface  times  diffusion  length.  Both  interface  surface  and  diffusion  length  are 
increasing  functions  of  time.  Having  these  two  functions  explicitly  one  can  calculate  a  reacting  volume  formation  and,  with  chemical  kinetics 
scheme  known,  the  combustion  process  can  be  described  in  detailed.  Such  an  approach  has  been  developed  in  IVTAN  for  simulation  of  high¬ 
speed  flow  of  chemically  reacting  mixture  for  plasma  chemistry  applications  and  simulation  of  combustion  in  compressor-type  chemical 
reactors  basing  on  Diesel  engine. 

Because  of  very  short  residence  time  is  typical  of  ram/scram  jet  operation  conditions  the  measures  providing  a  significant  increasing  rate  of 
reacting  volume  are  very  actual. 

The  proposed  Advanced  Mixing  approach  is  to  increase  the  reacting  volume  growth  by  (1)  pulse  electrical  discharge  -  JE  energy  release  and 
(2)  by  MHD  interaction  -  J?B  body  force. 

Furthermore  it  was  found  recently  both  experimentally  and  theoretically  [1]  that  under  certain  conditions  typical  of  so  called  Plasma 
Aerodynamics  experiments  the  filamentary  structure  of  high  frequency  discharge  (streamer  discharge)  results  some  times  in  high  speed 
counterflow  jet  formation.  Such  a  jet  can  transfer  a  significant  mass  of  fluid  from  one  point  to  another  rather  distant  location.  Thus,  providing  a 
controllable  filamentary  discharge  across  a  primary  fuel/oxidant  interface  one  can  intensify  the  initial  mixing  of  fuel/oxidant  composition. 

The  second  component  of  Advanced  Mixing  is  vortex  generation  by  a  highly  non-potential  J?B  body  force.  This  force  is  defined  by  current 
density  component  J?  perpendicular  to  the  external  applied  magnetic  field  B. 

The  preliminary  information  about  such  a  mixing  process  was  collected  in  IVTAN  while  an  MHD  generator  flow  with  current-currying  non¬ 
uniformities  was  studied.  At  that  time  the  main  question  was  how  to  decrease  the  mixing  process  intensified  by  vorticity  generation.  Now  the 
amplification  of  mixing  by  MHD  and  Pulse  Plasma  affects  should  be  studied  in  turn  and  recommended. 

We  have  developed  earlier  (  see,  for  example  [2,3]  and  others  publications)  physical  and  numerical  models  to  treat  the  kinematic  and  diffusion 
mixing  in  gas-plasma  flows  to  evaluate  the  life  time  of  plasmoids  in  such  a  flow  under  intensive  MHD  interaction  conditions. 

It  is  proposed  now  at  first  to  adopt  these  approaches  to  the  physically  similar  situation  of  mixing  of  gaseous  fuel  and  oxidant. 

The  second  goal  of  the  project  is  to  use  the  models  developed  under  the  first  task  to  conduct  parametric  study  of  effects  of  Advanced  Mixing 
on  characteristics  of  fuel/oxidant  mixing  in  high  speed  flows. 


The  third  goal  is  the  project  is  to  formulate  preliminary  the  scheme  of  the  experimental  validation  of  proposed  approach  under  experimental 
conditions  available  at  IHT  RAS. 
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I.  INTRODUCTION 


One  of  the  main  problems  of  stable  combustion  in  high-speed  flow  is  an  effective  mixing  of 
gaseous  fuel  and  oxidant  providing  full  combustion  during  the  resident  time  in  combustor  chamber. 
Generally  speaking  in-flow  mixing  consists  of  two  main  mechanisms:  kinematics  mixing 
(convection)  and  molecular  diffusion.  The  latter  is  a  final  stage  of  mixing  needed  to  provide 
extended  combustible  mixture.  Kinematics  mixing  is  more  effective  in  vortex  flows;  the  higher 
vorticity  the  better  kinematics  mixing  which  means  the  larger  fuel/oxidant  interface  surface.  The 
second  mixing  mechanisms  -  diffusion  forms  a  reacting  volume  estimated  as  the  interface  surface 
times  diffusion  length.  The  only  place  where  the  combustion  can  occur  is  the  reacting  volume.  Both 
interface  surface  and  diffusion  length  are  increasing  functions  of  time.  Having  these  two  functions 
explicitly  one  can  calculate  a  reacting  volume  formation  and,  with  chemical  kinetics  scheme 
known,  the  combustion  process  can  be  described  in  detail.  Such  an  approach  has  been  developed  in 
IVTAN  for  simulation  of  high-speed  flow  of  chemically  reacting  mixture  for  plasma  chemistry 
applications  and  simulation  of  combustion  in  compressor-type  chemical  reactors  basing  on  Diesel 
engine. 

Because  very  short  residence  time  is  typical  of  ram/scram  jet  operation  conditions  the 
measures  providing  a  significant  increasing  rate  of  reacting  volume  are  very  actual. 

The  proposed  Advanced  Mixing  approach  is  to  increase  the  reacting  volume  growth  by  (1) 
pulse  electrical  discharge  and  (2)  by  MHD  interaction  -  JxB  body  force. 

It  was  found  recently  both  experimentally  and  theoretically  that  under  certain  conditions 
typical  of  so  called  Plasma  Aerodynamics  experiments  the  filamentary  structure  of  high  frequency 
discharge  (streamer  discharge)  results  sometimes  in  high  speed  jet  formation.  Such  a  jet  can  transfer 
a  significant  mass  of  fluid  from  one  point  to  another  rather  distant  location.  Thus,  providing  a 
controllable  filamentary  discharge  across  a  primary  fuel/oxidant  interface  one  can  intensify  the 
initial  mixing  of  fuel/oxidant  composition. 

The  second  component  of  Advanced  Mixing  is  vortex  generation  by  a  highly  non-potential 
JxB  body  force.  This  force  is  defined  by  current  density  component  J_l  perpendicular  to  the  external 
applied  magnetic  field  B. 

The  preliminary  information  about  such  a  mixing  process  was  collected  in  IVTAN  while  an 
MHD  generator  flow  with  current-currying  non-uniformities  was  studied.  At  that  time  the  main 
question  was  how  to  decrease  the  mixing  process  intensified  by  vorticity  generation.  Now  the 
amplification  of  mixing  by  MHD  and  Pulse  Plasma  affects  should  be  studied  in  turn  and 
recommended. 
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In  practice  the  formation  of  reacting  volume  is  provided,  the  first  -  by  diffusion,  which  is 
rather  conservative  and  under  conditions  of  interest  slow  process,  and,  the  second  -  by  kinematics 
mixing.  The  latter  can  be  very  intensive  due  to  turbulence  generated  at  the  unstable  interface  of  free 
jet  in  a  resting  media.  However,  even  in  such  a  case  the  significant  length  is  needed  to  provide 
kinematics  pre-mixing  of  whole  fuel-oxidizer  mass  flow.  Furthermore,  the  kinematics  pre-mixing 
resulting  in  a  multi-scale  irregular  fuel-oxidizer  structure  consisting  of  stretched  and  deformed 
original  pure  fuel  and  oxidizer  volumes  is  not  yet  enough  to  provide  a  combustion  conditions.  Only 
diffusion  leads  to  a  combustible  mixture  formation.  The  role  of  kinematics  mixing  is  to  reduce  of 
characteristics  scale  (thickness)  of  inter-component  stratified  structure.  It  seems  that  the  kinematics 
mixing  stage  required  for  the  preparation  of  a  combustible  mixture  is  defined  by  the  condition  when 
the  stratified  system  characteristics  length  Is  is  comparable  with  the  diffusion  length.  The  diffusion 
length  time  evolution  can  reliably  be  estimated  with  a  square  root  law.  Thus,  the  diffusion  time 
required  for  the  desirable  mixing  is  estimated  as  to  ~  h  ^/D.  For  a  rum  jet  combustor  applications 
this  estimation  is  followed  by  the  condition  that  requires  the  mixing  time  to  be  less  than  the 
residence  time  tR  ~  Lod/c  with  Lc  and  Vc  stand  for  the  combustor  length  and  in-combustor  velocity, 
correspondingly.  Thus,  the  practical  situation  analysis  came  back  to  the  study  of  fundamental 
problem  proposed  for  this  project  -  the  Advanced  mixing  utilizing  non-mechanical  mechanisms  to 
intensify  the  reacting  volume  creation.  The  additional  problem  related  to  feature  is  to  check  the 
possible  effects  of  electric  field  on  the  kinetics  of  combustion  and,  in  particularly,  on  the  ignition. 
The  expected  effect  is  to  improve  (accelerate)  the  ignition  stage  due  to  non-thermal  formation  of 
active  radicals  by  the  strong  electric  field  applied  or  arisen  in  the  flow  through  MHD  interaction. 
Such  an  extra  thermal  formation  of  the  radicals  (e,  O,  OH,  OH,  and  others)  can  be  provided  by 
the  electric  field  directly  in  ionization  and  dissociation  induced  by  electron-molecule  collisions. 

It  is  clearly  seen  from  above  brief  description  of  the  process  considered  here  that  the  very 
sophisticated  physical  and  mathematical  models  are  needed  to  treat  these  phenomena  at  acceptable 
accuracy.  From  the  other  hand,  not  all  input  data  required  for  the  closed  formal  description  are 
available,  and  for  this  reason  the  experimental  validation  is  absolutely  necessary. 

Basing  on  the  preliminary  results  of  analysis  conducted  in  the  frame  of  the  preceding  work 
under  contract  with  FOARD  [1]  the  approach  of  this  study  consists  of:  (1)  -  the  mathematical 
model  development  utilizing  the  full  Navier-Stokes  equation,  reduced  Maxwell  equation  to  treat  the 
MHD  interaction  process,  and  the  chemical  kinetics  solver  with  a  appropriate  chemical  kinetics 
scheme;  (2)  -  parametric  study  with  the  numerical  model  developed;  and  (3)  -  the  experimental 
setup  design  development.  In  order  to  reduce  the  computational  time  defined  mainly  by  the 
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chemical  kinetics  scheme  used  in  this  study  a  rather  compact  scheme  of  hydrogen-oxygen 
combustion  is  applied. 

In  fact  the  main  theoretical  and  numerical  models  needed  have  been  developed  earlier. 

In  the  present  work  the  following  physical  effects  are  of  primary  interest: 

a)  Molecular  viscosity  and  heat  conductivity.  The  turbulence  effects  are  hoped  to  be  resolved 
directly  rather  then  via  some  turbulence  model.  Most  of  turbulence  models  are  capable  of 
describing  the  mean  flow  only  if  the  flow  is  stationary  or  the  flow  field  varies  slowly.  We 
shall  consider  the  earlier  stages  of  mixing  process,  when  undeveloped  turbulence  may  exist. 
At  this  stage  we  consider  originally  undisturbed  flow  stage,  therefore  the  only  mechanism  of 
mixing  seems  to  be  the  molecular  transport. 

b)  Multi-component  diffusion  Traditional  approach  to  the  diffusion  of  particles  in  multi- 
component  mixtures  will  be  applied  in  this  work.  It  is  assumed  that  the  diffusive  flux  of  any 
species  is  proportional  to  the  local  gradient  of  its  concentration  (Pick’s  law).  The  baro- 
diffusion  and  thermo-diffusion  fluxes  can  be  neglected  in  the  flows  under  consideration. 

c)  Magneto-hydrodynamic  interaction.  This  is  the  main  effect  driving  the  flow.  The  arc-wise 
plasma  formation  can  be  created  either  in  fuel  or  oxidizer  and  sustained  due  to  external 
electric  source.  The  arc  current  interacts  with  the  transversal  magnetic  field  so  that  the  arc 
moves  toward  the  fuel-oxidizer  interface  and  enhances  stirring  of  components.  For  the  goals 
of  the  work  the  quasi-steady  approach  to  MHD-interaction  is  certainly  applicable:  the 
characteristic  speeds  of  the  processes  under  consideration  are  much  smaller  compared  to 
speed  of  light.  In  addition,  the  induced  magnetic  field  due  to  electric  currents  can  be 
neglected,  so  the  current  flow  under  constant  external  magnetic  field  can  be  considered.  The 
approach  to  description  of  plasma  properties,  of  which  the  electrical  conductivity  is  of  primary 
importance,  is  as  follows.  We  expect  that  for  low-temperature  (several  thousand  degrees  of 
Kelvin)  poorly  ionized  plasma  the  approximations  based  on  the  Saha-equilibrium  model  are 
reasonable;  so  that  the  conductivity  can  be  considered  as  the  function  of  gas  temperature. 
Same,  in  principal,  could  be  valid  for  so-called  Hall  parameter.  The  radiation  from  the  plasma 
will  be  taken  into  account  via  approximation  of  radiative  power  as  function  of  temperature. 

d)  Finite-rate  chemistry.  For  the  reacting  flow  computations  we  shall  apply  a  general  finite-rate 
chemistry  approach,  based  on  the  solution  of  chemical  kinetics  equations  to  obtain  the 
species  rates  of  production  at  each  spatial  point.  Any  chemical  reaction  is  assumed  to 
subject  to  the  acting-mass  law  and  reaction  rate  constants  are  assumed  to  be  of  Arrhenius 
type. 
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In  this  report,  the  mathematical  model  including  the  complete  set  of  governing  equations  is 
presented  in  Section  II,  the  computational  model  and  solution  algorithm  are  described  in  Section  III. 
Specific  features  of  the  numerical  experiment  setup  are  discussed  in  Section  IV.  The  main  results  of 
parametric  study  are  presented  in  Section  V.  The  initial  stage  of  experimental  facility  development 
is  briefly  described  in  Section  V.  The  summary  of  this  study  is  presented  in  the  Concluding 
Remarks. 
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II.  MATHEMATICAL  MODEL 


The  flows  under  consideration  can  be  described  by  the  coupled  system  of  two-dimensional 
equations  reflecting  the  conservation  of  mass,  momentum  and  total  energy  for  the  whole  mixture 
(Navier-Stokes  equations),  conservation  of  mass  of  different  species  with  taking  in  to  account  the 
chemical  conversions.  The  system  is  enclosed  by  the  equations  of  state,  by  definition  of  transport 
coefficients,  by  definitions  of  chemical  rates  of  productions  for  every  species  and  by  relationships 
describing  the  electromagnetic  interaction.  Below,  the  general  mathematical  model  is  given  as  it  is 
implemented  in  the  code. 


au  ar  ac  ar^  ac^  ^ 

—  -I-  —  -I-  —  -I- — ^-1- — ^  =  s 
dt  dx  dy  dx  dy 


(1) 


U  is  the  conservative-variable  vector. 


f  P 


PVy 

pE 

Pn 


[pYn 


) 


(2) 


where  p  is  the  fluid  density,  Vx  and  Vy  are  the  velocity  components,  E  is  the  total  specific  energy,  and 

N 

{Yi,  i=l,N}  are  the  species  mass-fractions,  T=p//p,  52^  “ 

i  =  l 

The  inviscid  flux-vectors  are  given  as 


f  PVx  ] 

^  PVy  " 

P^x  +  P 

PVx^y 

P^x^y 

Pv^+P 

VxPE  +  V^P 

,G  = 

VypE+VyP 

pVxYl 

PVyYl 

y  P^xYn  J 
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Here,  P  thermodynamic  pressure  defined  later. 

Diffusive  fluxes  are  described  as  follows: 


(  0  ^ 

^  0 

xy 

xy 

'^yy 

^x 

’  G,,  - 

ey 

. . 

•f  y\ 

^JyN 

(4) 


The  components  of  viscous  stress  tensor: 


dv  2 


dx 


dv,. 


-+■ 


dy  dx 


2 


(5) 


dv  ^Vy 

where  p.  is  the  molecular  viscosity,  Vv  = — --\ - The  Stokes’  hypothesis  on  the  second 

dx  dy 

viscosity  has  been  taken  into  account  in  (5). 


Energy  fluxes  are  specified  as: 


-  ( 

^x~  '^Vx'^xx '^^y'^xyj  '  XI 

ax 


-  ^ 
j  ^  dy  ^y'^yy 


i=l 

)+Y,hiJ  yi 


(6) 


i=l 


Here,  A,  is  the  heat  conductivity  and  T  is  the  temperature.  The  last  terms  in  (6)  represent  the 
enthalpy  fluxes  due  to  multi-component  diffusion. 

Diffusive  fluxes  are  defined  according  to  Pick’s  law: 
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xi  P^im  ■\  ’  ^ 


yi  P^ii 


dx 


,  i=l,N 


(7) 


Transport  coefficients  Dim,  |J.,  A.. 


Di^^a-Xi) 


(8) 


Xi-W  ■  Yi jWi  is  the  mole  fraction  of  i-th  species,  W,  is  the  molecular  weight  of  i-th  species  and 


'  N 


W  =l/52^i  A^i  molecular  weight  of  mixture.  The  binary  diffusion  coefficients.  Ay,  are 

/  (=1 

defined  as  follows: 


D^.  =  2.68. 


10 


-7 


'^’^^{Wi+Wj)l{2wyvj) 


m 

s 


(9) 


1  H= 

Here,  =  — (a,  -i-a^  )  represent  the  characteristic  collision  diameter,  Tij  -kTjeij  (k  is  Boltzman 


constant)  is  the  characteristic  temperature,  8^  =  .^8,  •  8^  is  the  parameter  of  function  describing  the 
inter-molecular  interaction. 

Collision  integral  is  approximated  as  =1.074-(7)y 

Viscous  coefficient,  p,,  is  calculated  as 


N 

(=1 


7 

j*i  j 


where  /-th  species  viscous  coefficient  is  defined  as 


(10) 
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(11) 


(1/ 


2.6693  10  ^ 


kg 
ru' s 


Collision  integral  np'^^(7]*)  is  approximated  as  =  1.157  •  (T*)  ,  and  functions  G,)  are 

defined  as 


The  heat  conductance  coefficient  is  defined  as  follows: 


N 


X=ZK 

i=\ 


1  +  1.0651:  G^.-^ 

V  j*'  '  y 


X-l 


where  X-  =  X,-  (o.  1 15  +  0.354  •  C^,-  //?, ) 


(12) 


functions  G.-,-  are  defined  as 


[!  +  (>.,/>., TV, 


/-th  species  heat  conductivity  is  defines  as 


r\ 


1  ^ 

I'W, 


(13) 


is  the  universal  gas  constant,  /?,  is  the  gas  constant  of  i-th  species,  Ri=I^IWi,  Cn  and  Cpi  are  species 
heats  at  constant  volume  and  constant  pressure,  respectively. 
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Thermodynamic  properties  of  mixture. 


Equation  of  state 


P  =  pRT 


(14) 


where  P  is  the  pressure,  T  is  the  temperature,  R  is  the  gas  constant,  R=I^/W. 


Caloric  relation 


e  =  h-  P/p 


(15) 


where  e  is  the  specific  internal  energy  and  h  is  the  specific  enthalpy.  Total  energy  is  defined  as 


E-e  +  -{vl+Vy)  (16) 

and 

h^^YMT)  (17) 

(=1 


i-th  species  enthalpy,  hi(T),  is  expressed  as 


h{T)^k,f  + 


T 

icpi(T)dT 


(18) 


where  hij  is  the  formation  enthalpy  of  i-th  component. 


The  source-term  vector  S  is  given  as 
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Here,  Sx  and  Sy  represent  the  components  of  electro-magnetic  force,  Se  represents  the  energy  rate  of 
change  due  to  Joule  heating  and  radiation.  These  terms  are  described  in  the  sub-section  “Plasma 
Model”  below.  The  terms  {w,  },  i=l,N  represent  the  rates  of  production  of  species  due  to  chemical 
reactions.  They  are  described  in  “Finite-Rate  Chemistry  Model”. 

Plasma  Model. 

In  general,  the  electrodynamics  of  steady-state  weekly-ionized  gas  discharge  can  be  described  by 
Ohm’s  law: 

j  +  +  (20) 

B 

where  j  is  the  electric  current  density,  E  is  the  electric  field  strength,  a  is  the  electric  conductivity, 
B  is  the  magnetic  field  induction,  and  p  is  the  Hall  parameter. 

The  components  of  electro-magnetic  force  are  defined  as 

S,,,  =  [JxBL,,  (21) 

Energy  rate  of  change  is 
“  (j  ’E)-  Qrad 

where  approximation  to  the  radiative  power  Qrad  is  taken  in  the  form: 
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(23) 


Qrad  - 


where  £  is  the  effective  absorption  coefficient,  and  Osb  is  the  Stefan- Boltzman  constant. 

Coupled  with  the  definition  of  electric  potential,  E=-grad  tp,  the  Equ(20)  can  be  reduced  to 
single  two-dimensional  (three-dimensional  in  general)  equation  for  potential. 

Finite-Rate  Chemistry  Model 

The  chemical  kinetics  model  is  applied  to  determine  the  species  rates  of  production  {w,- }  due 
to  chemical  conversions,  or  reactions. 

Any  chemical  reaction  can  be  represented  symbolically  as: 


(24) 


where  v'lr  and  V'lr  are  the  stoichiometric  coefficients  of  reagents  and  products,  respectively. 

The  i-th  species  rate  of  production  is  then  defined  in  accordance  with  the  acting-mass  law: 


Here,  ci=p-Yi/Wi  are  the  molar  concentrations,  kfr  and  kbr  are  the  reaction  rate  constants,  forward  and 
backward,  respectively.  The  Arrhenius  form  is  assumed  for  both  kjr  and  kbr- 

kf^^a-T'’  ■  exp (-  E^/T)  (26) 

where  a,b  and  Ea  are  reaction  parameters.  These  can  be  evaluated  theoretically  or  from 
experimental  data.  For  the  hydrogen-air  mixture  these  data  are  tabulated  as  well  as  that  for  many 
other  reactions. 
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The  Model  Up-grade 


The  first  stage  study  with  the  model  described  above  has  concerned  the  feasibility  of  the  electric  arc 
to  initiate  and  maintain  the  combustion  via  the  motion  across  the  fuel/oxidizer  interface  due  to 
electromagnetic  field.  It  was  numerically  shown  that  such  a  discharge  both  initiates  combustion 
process  in  the  mixing  layer  due  to  high  temperature  of  the  arc  and  increases  the  mixing/combustion 
volume  due  to  arc  motion  and  enhanced  diffusion  of  species  in  the  volume.  The  combustion  was 
detected  to  be  self-sustained  during  rather  long  time  even  the  current  and/or  magnetic  field  is 
switched  off. 

In  those  computations  rather  simple  models  of  the  medium  were  used. 

First,  the  chemical  kinetics  mechanism  used  is  considered  to  be  appropriate  to  the 
combustion  under  thermal  equilibrium.  In  addition,  the  kinetics  under  temperatures  higher  than 
3000K  can  differ  from  that  considered  before.  Therefore  sophisticated  kinetic  scheme  should  be 
tested  which  includes  interactions  with  the  charged  particles,  first  of  all,  electrons.  The  neutral- 
neutral  interactions  should  be  revised  for  temperatures  higher  than  3000K. 

Second,  the  transport  properties  of  the  multi-component  mixture  should  be  revised  due  to 
availability  of  electric  field  and  charged  particles.  First  of  all  the  conductivity  of  plasma  should  be 
computed  more  accurately.  The  electron-neutral  and  electron-ion  collisions  are  of  importance.  In 
general  these  interactions  depend  upon  electron  concentration  and  electron  temperature,  the  latter 
being  different  from  heavy  particles  temperature.  Other  transport  properties  such  as  viscosity,  heat 
conductivity  and  diffusion  coefficients  can  significantly  be  effected  by  the  charged  particles 
interactions.  Moreover,  the  approach  to  mass-diffusion  used  up  so  far  should  be  re-evaluated.  This 
approximate  approach  is  based  on  the  Tick’s  law  and  effective  diffusion  coefficients.  The  approach 
is  known  to  be  stable  while  time-integration,  but  may  be  inadequate  in  plasmas  with  the  deep 
property  gradients.  In  this  case  the  solution  Stefan-Maxwell  equations  may  become  necessary. 

Third,  the  radiation  emission  from  the  discharge  region  was  estimated  with  the  approximate 
model.  It  was  assumed  that  total  radiation  losses  from  optically  thin  plasma  can  be  calculated  by 
simple  law  taking  into  account  the  local  absorption  coefficient.  Special  study  is  still  needed  to 
understand  whether  the  radiation  losses  are  important  from  the  viewpoint  of  energy  balance.  If  so, 
they  should  be  evaluated  with  more  sophisticated  model. 

The  following  sections  are  primarily  intended  to  analyze  the  plasma  models  used  before  and 
to  select  the  models  more  suitable  for  the  discharges  under  consideration.  Three  items  will  be 
analyzed  below:  chemical  kinetics,  transport  properties  and  radiation  losses. 

In  Table  1  a  typical  chemical  kinetics  scheme  is  presented.  This  scheme  was  used  in  the 


first  series  of  numerical  simulation. 
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Table  1.  List  of  reactions  for  H2-O2  combustion 


1)  H2+02^0H  +  OH 

2)  OH  +H2^H20  +  H 

3)  H  +  02^0H  +  0 

4)  0  +  H2^0H  +H 

5)  O  +  H2O  <^OH  +OH 

6)  H  +  H  +  M  ^H2+M 

7)  H  +  O2+M  ^H02  +  M 

8)  OH  +H  +M  ^H20  +  M 

9)  HO2  +  H2<^  H2O2  +  77 

10)  HO2  +  7/O2  <^H202+  O2 

11)  H  +  H02^0H  +  OH 

12)  H  +  H02^H20  +  0 

13)  H  +  H02^H2+02 

14)  0  +  H02^0H  +  O2 

15)  OH  +H02<^H20  +  02 

16)  OH  +OH  +  M  ^H202+M 

17)  HO2  +  H2<^  H2O  +  OH 

1 8)  HO2  +  H2O  H2O2  +  OH 

19)  H  +  H202  ^H202  +  0H 

20)  OH  +M  ^0  +  H  +  M 

21)  0  +  0  +  M<^02+M 

It  takes  into  account  all  main  stages  of  the  combustion  process:  excitation,  ignition  and  heat 
release  stage.  The  mixture  is  assumed  to  be  in  local  thermal  equilibrium. 

However,  under  plasma  conditions,  the  electron-impact  reactions  may  play  the  significant 
role.  For  this  reason  the  first  simple  chemical  kinetics  scheme  of  Table  1  is  now  reconsidered  to  be 
fit  to  the  conditions  of  the  experiments  planned  for  future  study. 

The  following  interactions  are  planned  to  be  included  in  the  scheme. 
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1)  Dissociation  reactions 

H2  //+//  +  e 

O2  +  s  0  +  0  +  6? 

OH  +e^O  +  H  +e 

2)  Major  ionization  reactions 
H2  ^  H2  +6  +  6 

O2  +  6  O^  +6  +  6 
H  +  e  +  e  +  e 

0  +  e^  0'*'  +  6  +  6 
OH  +  6  OH "''  +  6  +  6 

This  set  of  reactions  is  thought  to  be  important  on  any  stage  of  the  combustion  process.  The 
data  on  these  reactions  are  now  being  sought  in  the  literature.  The  role  of  direct 
dissociation/ionization  of  water  is  now  investigated.  Interactions  of  electrons  with  several 
intermediate  species  like  HO2,  H2O2  will  probably  be  neglected  since  their  life-time  is  rather  small 
as  well  as  their  concentrations. 

Ionization  of  main  species  due  to  neutral-neutral  collisions  of  type  of 

X+M^X++M+6 

where  denotes  It,  O2,  O,  H,  OH,  and  M  denotes  any  of  the  species  will  probably  be  not  considered 
since  such  reactions  occur  under  high  temperatures  (probably  higher  that  leV).  At  the  same  time 
electron-impact  reactions  may  take  place  due  to  high  electron  temperatures,  which,  in  turn,  may 
take  place  in  high  electric -field-strength  discharge.  Therefore,  the  set  of  dissociation/ionization 
reactions  formulated  above  should  depend  on  the  electron  temperature. 

Another  aspect  of  the  scheme  being  developed  arises  from  consideration  of  higher- 
temperature  range.  Within  the  arc,  the  temperatures  higher  than  3000K  may  occur.  Hence,  the 
dissociation  reactions  of  type  of 

H2  +  M  <^H  +  H  +  M 

which  were  included  in  the  previous  scheme  may  occur  with  the  reaction  rates  different  from  those 
considered  before.  The  preliminary  list  of  reactions  which  are  assumed  to  be  included  in  the  kinetic 
scheme  is  given  in  Table  2. 
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Table  2.  Suggested  list  of  reactions  for  H2-O2  combustion  with  ionization 


1)  H2+e^H+H  +  e 

2)  O2  +  s O -\-0 c 

3)  OH  +  6  O  H  c 

4)  H  +  e  +  e  +  e 

5)  O  +  e  O^  +  e  +  e 

6)  OH  +e<^OH^  +e  +  e 

7)  H2~^  C  H2  S -h  C 

8)  O2  +  ^  O2  +  £?  +  ^ 

9)  H2+02^0H  +OH 

10)  OH  +H2^H20+H 

11)  H  +  02^0H  +0 

12)  0+H2^0H  +  H 

13)  0  +  H20^0H  +  OH 

14)  H  +  H+M^H2  +  M 

15)  H  +  O2+M  ^H02+M 

16)  OH  +H  +  M  ^H20  +  M 

17)  H02+H2^H202+H 

1 8)  //O2  +  HO 2  ^  H2O2  +  O2 

19)  H  +  H02^0H  +OH 

20)  H  +  H02^H20  +  0 

21)  H  +  H02^H2+02 

22)  0+H02^0H  +O2 

23)  OH  +  HO2  <^H20  +  02 

24)  OH +OH +M  ^H202+M 

25)  HO2  +  H2^  H2O  +  OH 

26)  HO 2  +  H2O  H2O2  +  077 

27)  //  +  H2O2  /72O2  +  OH 

28)  0//+M^0  +  //+M 

29)  0  +  0  +  M<^02  +  M 

Transport  properties 


In  the  previous  work  the  transfer  of  species,  viscous  momentum  and  heat  was  modeled  by  the  well- 
known  multi -component  diffusion  model  formulated  in  [2].  Commonly  used  simplification  of 
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multi-component  diffusion  applied  the  effective  diffusion  coefficients  and  Pick’s  law  to  estimate 
species  mass-flux  was  implemented  in  the  code.  All  transport  coefficients  were  estimated  with 
using  so  called  eollision-integrals  of  viseous  and  diffusion  types.  In  turn,  the  collision-integrals  are 
obtained  by  integration  of  collision  cross-sections  over  the  energy  space.  The  values  of  collision- 
integrals,  fi(/,s),  were  taken  as  approximations  of  temperature  functions  [3].  These  approximations 
used  the  Lennard-Jones  interaction  potential  to  obtain  the  deflection  angle  of  the  eolliding  particles 
as  function  of  impact  parameter  and  relation  energy,  then  the  relevant  cross  section,  and,  finally,  the 
collision  integrals.  The  Lennard-Jones  interaction  potential  model  is  known  to  work  well  for 
temperatures  below  nearly  lOOOK.  At  higher  temperatures  another  models  should  be  used  for 
interaction  potential.  In  addition,  the  interaction  of  charged  particles  should  take  into  account  the 
sereened  Coulomb  potential. 

Another  aspect  of  the  problem  is  related  to  using  the  Pick’s  law  for  mass-flux  expression.  It 
is  easy  to  implement  into  the  eode,  but  the  accuracy  of  the  approach  remains  under  question.  More 
relevant,  although  time-consuming,  is  the  solution  of  Stefan-Maxwell  relations  for  mass-diffusion 
fluxes.  The  set  of  Stefan-Maxwell  equations  ean  be  expressed  as 


XjJj 

MjD^j 


MAj 


,i=l,N-l 


(27) 


where  7,  is  the  mass-flux  of  species  i  due  to  diffusion  (only  concentration  diffusion  is  taken  into 
account),  x,  is  the  mole  fraetion  of  species  i,  M,  is  the  moleeular  weight  of  i-th  species  and  A>  are 
the  binary  diffusion  coefficient.  System  (27)  should  be  resolved  for  7.  This  requires  matrix 

inverting,  J  =  D“^(Vx),  which  is  extremely  costly.  Therefore,  iterations  are  usually  applied  to 

obtain  7.  Nevertheless,  this  procedure  still  remains  expensive  compared  with  simplified  models. 
Probably  one  of  the  best  from  the  viewpoint  of  cost-accuracy  consideration  is  the  Yos  model.  The 
model  gives  diffusion  flux  of  species  i  as 


7  Mi 

*  ^  M 


1-x, 


V 


£).  -  ^  ^ 
^im  — 


(28) 
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Here  Am  is  the  effective  diffusion  coefficient,  c,  is  the  mass-fraction  of  species  i. 

In  the  case  of  Stefan-Maxwell  relations  the  driving  force  can  take  electric  force  into  account 
directly.  In  Yos  model  it  is  neglected.  Regardless  of  the  diffusion  model,  the  binary  diffusion 
coefficients  have  to  be  calculated.  Their  calculation  requires  the  knowledge  of  collision  integrals. 
Two  approaches  are  planned  to  try  for  collision-integral  calculations. 

First,  collision-integral  for  neutral-neutral  collisions  is  approximated  as  function  of 

temperature:  =  [a^j  -l-^y  Inj)^,  where  parameters  a,)  and  bij  can  be  found  from  known  values 

of  (T1=300K)  and  (T2=20000K)  and  assuming  that  Lennard-Jones  interaction 

potential  is  valid  at  low  temperatures. 
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^  =  0A05rl]n{\  +  2A)  [A^], 


if  1<A<100 


and 


=  ^+- ^  =  O-S^L  (in  A  -  0.305) , 


if  A>100. 


Here 


rl{T)- 


(  1  f 

AtI£^T 


[A^], 


ro  .T 


312 


A(r,nA  =  ^  =  A 


1/2  • 


To  caleulate  viseosity  and  heat  eonduetivity  we  plan  to  use  approximations  of  type  of 
Wilke- Vasil’eva  [7]: 


N 


5Cj- 
i=\ 


^  A/f 


M: 


c„:  +2.5/?®  (1.55c,  -1)— 


where 


rf; 


The  Sehmidt  number  5c, (7)  =  — .  It  is  assumed  that  for  all  eomponents. 

■ 

Some  data  were  also  taken  from  [8]. 
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In  second  approach  the  collision-integrals  values  can  be  taken  in  closed  form  as  functions  of 
temperature  [9].  In  this  work  the  collision-integrals  were  calculated  for  many  equilibrium  air 
species  interaction  including  the  charged  particles  interactions. 

In  the  future  it  is  assumed  to  implement  both  transport  coefficients  models  to  take  into 
account  charged  particles  interactions. 

Radiation  from  plasma 

In  the  calculations  made  earlier  the  simplest  radiation  model  was  used.  The  model  used  the 
concept  of  “gray-body”  with  effective  absorption  coefficient.  Such  concept  just  estimates  the  level 
of  radiation  power,  probably  upper  limit  of  it.  Currently,  we  don’t  know  how  accurate  are  these 
estimations.  Detailed  analysis  of  radiation  emission  is  extremely  complex  problem,  and  it  is  hardly 
possible  to  do  by  ourselves.  Nevertheless,  work  is  now  carried  out  to  find  acceptable  level  of 
radiative  losses  description. 

Concluding  remarks  on  the  mathematical  model  development 

The  time  period  of  last  several  months  has  been  devoted  to  the  analysis  of  the  models 
applied  to  the  investigations  of  MHD  assisted  mixing/combustion  in  non-premixed  fuel/oxidizer 
system.  From  examination  of  the  literature  we  drew  a  conclusion  that  some  modifications  and 
improvements  should  be  made  both  in  physical  and  computational  models  of  the  processes 
analyzed. 

First,  the  kinetic  mechanism  should  be  corrected  to  take  into  account  reactions  with  charged 
particles  and  to  expand  the  model  to  higher  temperature  range.  The  preliminary  list  of  reactions  has 
been  compiled,  and  now  the  work  is  done  to  implement  it  to  the  code  and  to  test  it. 

Second,  the  transport  properties  models  were  analyzed  since  we  think  that  the  current  model 
is  not  completely  adequate  to  simulate  the  transfer  in  non-equilibrium  system  with  the  charged 
particles  and  strong  property  gradients.  The  most  robust  candidate  seems  to  be  the  Stefan-Maxwell 
relations  for  mass-diffusion  fluxes.  At  the  same  time  this  model  is  known  to  be  rather  costly  from 
computational  point  of  view.  The  compromise  between  cost  and  accuracy  may  be  the  Yos  diffusion 
model.  Both  these  models  are  now  being  implemented  to  the  code  and  tested.  Also  modifications  to 
the  models  of  viscosity  and  heat  conductivity  are  being  introduced. 

Third  point  of  analysis  carried  out  seems  to  be  not  so  important  compared  to  first  two. 
Nevertheless  some  more  sophisticated  estimations  of  radiation  losses  from  the  arc  are  being  worked 
out. 
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m.  NUMERICAL  METHOD 


In  order  to  obtain  all  flow  fields  in  any  point  of  the  domain  under  consideration,  the  entire 
computational  domain  is  divided  onto  small  domains  called  cells.  Within  each  computational  cell, 
the  flow  variables  may  be  assumed  to  be  constant  (first  -  order  accurate  solution)  or,  for  example, 
to  be  linear  (second  -  order  accurate  solution).  Other  cell-distributions  are  possible,  but  are 
extremely  expensive  from  the  viewpoint  of  computational  time.  Only  constant  and  linear  variable 
distributions  are  implemented  in  the  code. 

The  system  (1)  is  discretized  in  space  and  time  starting  from  integral  conservation  laws 
(weak  formulation): 

-^jUdto+ JF^  +  (29) 

CO  aco  da 

This  is  integral  form  of  (1)  valid  for  any  computational  cell  with  the  volume  co  enclosed  by  the 
surface  Sa  (in  two-dimensional  case  to  is  the  surface  and  Sa  is  the  curve  surrounding  the  cell).  For 
both  constant  and  linear  variable  distribution-in-cell,  the  transient  term  in  (29)  can  be  approximated 
as  follows: 

—  \U  dco  =  co - , 

dt  J  At 

CO 

where  AC/=C/(t-i-At)-C/(t)  with  the  1-st  order  approximation  of  time-derivate.  The  surface  fluxes  are 
integrated  as  follows.  For  example, 

JF  =  Jf  +  Jg  dSl  =  j {Fn^  +  Gny)dS^^  =  (^  +  G-ny)AS^^  (30) 

Here,  AA®  is  the  surface  area  value,  Ux  and  Uy  are  the  components  of  unit  normal-vector  directed 
outward  the  cell  considered.  The  average-flux  value  depends  on  the  grid  partition  and  on  the 
variable  distribution-in-cell.  We  shall  consider  the  quadrilateral  cells,  so  that  any  internal 
computational  cell  has  four  neighbors.  The  relationships  (30)  are  then  rewritten  as  follows: 
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ye  CO  j 

where  Fj  stands  for  the  face-averaged  flux  and  Sj  stands  for  the  face  area.  The  face  implies  the  line- 
segment,  part  of  the  whole  cell-surface,  straddling  two  adjusting  cells.  Four  faces  enclose  each 
computational  cell.  For  F*-  and  2"'*  order  accurate  solutions,  the  one-integration-point  is  sufficient, 
i.e.  the  surface  flux  can  be  evaluated  at  the  mid-point  of  the  face: 

\FjdSj  =  F(U„)AS^.  =  [FiVJn,  +G(U„)njAS^.  (31) 

The  flux  components,  F  and  G,  were  presented  in  the  previous  Section.  The  dissipative  fluxes,  Fy 
and  Gy,  are  handled  with  analogously. 

The  source-term  is  discretized  as  follows 


jS  c/co  =  co  S(U(,3) . 


Hereafter,  all  the  flow  variables  imply  cell-volume-averaged, 

Jf/  c/®  =  CO-  U(Q,  or  simply,  ®-U  . 

Collecting  surface-flux  integrations  and  giving  variable  distribution-in-cell  leads  to  a  system  of 
algebraic  non-linear  equations  for  each  cell  state-vector  U: 

®— +  R  =  0  (32) 

At 


where  residual-vector  R  is  the  approximation  of  surface-flux  quadratures  plus  source  term: 


R 


[f  +  fS  da 


yeco 


(33) 


The  system  (32)  is  integrated  in  time  using  explicit  low-memory  Runge-Kutta  method: 


Uo  =  U" 
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(34) 


U,.=Uo-a,.— R(U,_i)  ,  /=!,/ 

CO 

=  U; 

where  U"  is  the  solution  on  time-level  f  and  is  the  solution  on  time-level 


The  time-step,  At,  is  selected  as  minimum  of  convective  and  diffusive  time-scales  for  all  the 
mesh  control-volumes.  While  source-terms  for  momentum  and  energy  equations  are  computed 
explicitly,  calculation  of  chemical  species  production  rates  is  carried  out  implicitly  employing 
special  chemical-kinetics  solver. 


Wi  = 


dipYj) 

dt 


chem 


]_ 

At 


f  +^t 

f  w,. 

t" 


dt 


(35) 


Integration  in  time  is  performed  by  the  program  KENEL  developed  in  IVTAN  Low- 
Temperature  Plasma  Division  for  solution  of  systems  of  ordinary  differential  equations  of  high 
stiffness  inherent  to  problems  with  chemical  and  thermal  non-equilibrium.  On  every  time-step.  At, 
system  (4c)  is  integrated  based  on  “chemical”  time-steps  and  implicit  treatment  of  production-rates. 
The  chemical  time-steps  are  evaluated  within  KENEL  on  the  base  of  eigen-values  of  chemical 
source-term  Jacobian  matrix.  Low-cost  LU-decomposition  of  Jacobian  is  applied  on  each  chemical 
time-step  to  solve  the  system. 

Implicit  treatment  of  chemical  source-term  allows  one  to  avoid  strong  restrictions  to  the  gas- 
dynamic  time-steps  imposed  by  finite-rate  chemistry. 

The  integration  of  inviscid  and  viscous  flux-vector  over  the  control-volume  surface  is 
performed  as  follows.  One-integration-point  on  the  each  surface-element  (cell-interface  straddling 
two  neighboring  cells)  is  suitable  to  obtain  both  1-st  order  solution  (piece-wise  constant  variable 
distribution  within  the  cell)  and  2-nd  order  solution  (piece-wise  linear  variable  distribution)  in 
space.  The  inviscid  flux  at  the  integration-point  is  found  from  exact  solution  of  Riemann  problem 
(Godunov’  method),  i.e. 


F,-p  —  F(Ujp) ,  and  - R/e(U^,U^) 


(36) 
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Rie^yJi,  Ur)  implies  solution  of  Riemann  problem  between  two  states  Ul  and  Ur,  left  and  right 
state-vectors  on  both  sides  of  the  interface.  For  two  adjacent  cells,  P  and  E, 

Ul  =  Up  -i-\|/p(VUp,Arp^) 

U„=U^->|«j{VU^,ArB,.)  (37) 

VUr,£  are  the  cell-centered  gradients  and  Alrm  is  the  vector  directed  from  cell-center  P{E)  to  mid¬ 
point  of  the  interface  M.  The  gradient  in  each  cell,  VUp,  is  calculated  by  the  least-square  method  as 
minimum  of  functional 

4.p=  £[A(/p^-(vUf.Arf,)f  , 

Je  cOp 

^Upj  =Uy-Up 

Here  J  stands  for  any  cell  from  the  stencil  associated  with  cell  P.  Usually,  five-point  stencil 
of  the  nearest  neighbors  (including  P)  is  considered.  The  limiter-functions  \|/p’s  are  taken  to  satisfy 
the  monotonical  distributions  near  the  strong  flow  discontinuities.  In  the  smooth-flow  regions,  the 
second-order  accurate  spatial  distributions  are  recovered.  In  one-dimensional  problem,  this  method 
becomes  equivalent  to  the  well-known  Van  Leer  MUSCL  method. 

Viscous  fluxes  at  cell  interface  may  be  constructed  in  two  ways.  In  first  way,  the  face- 
gradient  is  taken  to  be  the  algebraic  average  of  two  cell-gradients  described  just  above  with  adding 
up  the  projection  of  normal-derivative  to  prevent  the  non-physical  oscillations.  In  second  way,  the 
gradient  at  interface  is  calculated  by  application  of  Gauss’  theorem  to  auxiliary  cell 

®a  (^U^ )  =  J  (u  •  )  dS^  . 

3co^ 


Auxiliary  cell  represents  the  quadrilateral  which  vortices  are  two  cell-centers  {P  and  E)  and 
start-point  and  end-point  of  the  interface  segment  itself.  Both  approaches  were  found  to  be  equally 
appropriate. 
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IV.  PARAMETRIC  STUDY 


IV-1.  Arc-Driven  Mixing  and  Combustion 

To  highlight  the  main  effects  discussed  in  Introduction  we  consider  a  simplified  problem. 
Assume  that  the  rectangular  parallelepiped  is  filled  with  separated  air  and  hydrogen.  Let  the 
computational  plane  be  the  XT-plane.  The  size  of  the  domain  in  X-direction  is  taken  to  be  10  mm. 
The  size  of  the  domain  in  T-direction  is  also  taken  to  be  10  mm.  We  assume  that  variations  of  the 
flow  parameters  in  Z-direction  can  be  neglected,  so  that  the  two-dimensional  problem  in  XT-plane 
could  be  considered.  The  domain  is  filled  with  the  separated  fuel  components  in  such  a  way  that 
upper  half  of  the  domain  (T>5)  is  filled  with  air  and  lower  half  is  filled  with  hydrogen.  This  will 
referred  to  as  the  Case  1.  Alternate  filling  case,  i.e.  hydrogen  if  T>5  and  air  if  T<5,  will  be  referred 
to  as  the  Case  2.  At  the  initial  moment  the  arc  of  2  mm  in  diameter  is  ignited  so  that  the  center-point 
of  the  arc  locates  at  the  point  x=5mm,  y=3mm.  The  arc  column  is  specified  to  be  co-linear  to  Z- 
direction  and  is  assumed  to  maintain  from  the  external  electric  source.  In  all  the  calculations  the 
total  current  through  the  arc  10  Ampere  is  assumed  to  hold.  The  constant  magnetic  field  is  specified 
as  B  =  (5,0,0),  5  =  2  tesla.  In  such  formulation  the  vectors  of  current  density  and  electric  field  take 
the  following  form: 

J  =  (0,0,/,),  E  =  (0,0,5,). 

The  component  of  electric  field  vector  5,  is  the  same  everywhere  in  the  domain.  For 
simplicity,  we  neglect  with  the  induced  electric  field,  [VxB].  This  can  be  done  since  the 
characteristic  flow  speeds  of  the  order  of  100  m/s  were  expected  (and  calculated),  whereas  the 
applied  electric  field  varied  in  the  range  from  1  to  4  kV/m.  Then,  one  can  write  for  total  current 

^arc  ~  z  ~  \  ^  '  ~  ^z  '  ^^xy 

^Z  ~  ^arc/ ^^xy 

/,  =a(x,y)-5. 

Next,  the  electro-magnetic  force  has  only  one  component,  namely 
Sx  —  J  z '  B ,  Sy  =  0 
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The  energy  rate  of  change,  Se,  is  given  by  expression: 


S,=(JE)=7,.£, 

The  following  initial  conditions  are  specified.  Pressure  Pq  =  105  Pa,  temperature  Tq  =  300K  outside 
the  arc  and  Tq  =  6000K  inside  the  arc.  The  boundaries  of  the  domain  are  considered  to  be  the 
adiabatic  walls. 

Some  results  can  be  seen  from  Figures  1-3.  In  Fig.l,  the  no-chemistry  Case  1  is  presented. 
The  distributions  of  mass  fraction  of  hydrogen  (left)  and  temperature  (right)  are  shown  for  8  time- 
moments.  For  the  same  time-moments,  the  same  distributions  are  given  in  Fig.2  for  Case  1  (with 
chemistry).  Also,  the  main  combustion  product,  H2O,  is  shown  by  colored  contours.  Finally,  Case  2 
is  represented  in  Fig.3. 

The  first  observation,  which  immediately  can  be  done,  is  that  the  region  of  disturbance  of 
the  original  contact  boundary  is  essentially  higher  for  Case  1  compared  to  original  arc  diameter.  The 
characteristic  size  of  combustion  zone  is  comparable  with  the  size  of  the  domain.  The  combustion 
zone  is  well  detected  by  the  water  mass  fraction  distributions.  The  water  mass  fraction  reaches  the 
value  of  ~  0.21  (maximum)  at  the  final  stage.  The  maximum  temperature  varies  from  ~  5000K  at 
5mcs  to  3900K  at  the  final  moment.  Comparison  of  Fig.l  (no  chemistry  Case  1)  with  Fig.2  shows 
that  mixing  process  is  mainly  due  to  arc  motion.  In  both  cases  the  electro-magnetic  force  is  nearly 
the  same,  which  determines  the  similar  flow  patterns.  Although,  some  small  differences  exist:  the 
disturbance  zone  is  slightly  larger  for  the  reactive  case,  while  the  maximum  temperature  is  in 
average  higher  for  the  non-reactive  case. 

The  situation  becomes  quite  different  when  arc  is  ignited  in  the  air.  In  this  case  no  essential 
mixing  zone  enhancement  is  observed,  while  temperature  becomes  as  high  as  ~  7200K.  This  is  the 
fact  that  in  the  Case  1  the  arc  ignited  in  the  light  hydrogen  and  therefore  moving  quickly  “strikes” 
with  the  heavy  air,  which  acts  like  the  wall.  The  flow  reflects  from  the  contact  boundary  rather  then 
passes  through  it.  This  leads  to  smearing  of  the  arc  along  the  boundary  and  high  diffusion  (due  to 
high  temperatures)  provides  an  initiation  of  combustion  on  a  wide  front.  In  opposite  Case  2,  slowly 
moving  arc  ignited  in  the  heavy  air  begins  to  accelerate  when  coming  into  the  light  hydrogen. 
Therefore,  the  mixing  zone  surrounds  the  arc,  which  is  as  narrow  as  original  one.  In  the  Case  2,  the 
maximum  water  concentration  reaches  the  value  of  0.1.  It  is  clear  that  total  water  production  is 
significantly  higher  in  the  Case  1 .  Even  if  one  compares  the  flow  patterns  at  the  same  time  moments 
(35mcs),  the  more  developed  combustion  zone  in  Case  1  is  evident. 
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Fig.l.  The  mass  fractions  of  hydrogen  and  water(left)  and  temperature(right)  for  Case  1  (no 

chemistry). 
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Fig.la.  The  mass  fractions  of  hydrogen  and  water(left)  and  temperature(right)  for  Case  1  (no 

chemistry,  cont.). 
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Fig.2.  The  mass  fractions  of  hydrogen  and  water(left)  and  temperature(right)  for  Case  1. 
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Fig.2a.  The  mass  fractions  of  hydrogen  and  water(left)  and  temperature(right)  for  Case  1  (cont.). 


Final  Report 


32 


Final  Report 


33 


Fig.3a.  The  mass  fractions  of  hydrogen  and  water(left)  and  temperature(right)  for  Case  2  (cont.). 
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IV-2.  Post-Discharge  Reacting  Volume  Evolution 

In  the  previous  Section  the  numerical  study  of  the  enhancement  of  mixing/combustion  in 
hydrogen-air  mixture  by  the  MHD-driven  arc  column  has  been  carried  out.  In  these  experiments 
both  electric  current  and  magnetic  field  were  permanently  maintained  during  the  computations. 

In  order  to  understand  what  is  an  ignition  capability  of  such  an  MHD  driving  mechanism  the 
next  step  study  has  been  undertaken  to  observe  the  evolution  of  the  reacting  volume  after  the 
current  and/or  magnetic  field  switches  off.  The  problem  configuration  is  basically  the  same  as  it 
was  used  in  the  previous  Section.  The  following  two  problems  were  considered. 

Case  1:  as  well  as  in  the  previous  work  the  arc  initiates  in  hydrogen  and  enters  into  the  air  due  to 
electromagnetic  force.  25  microseconds  after  the  arc  generation  the  magnetic  field  is 
switched  off,  while  the  electric  current  remains. 

Case  2:  25  microseconds  after  the  arc  generation  both  electric  current  and  magnetic  field  are 
switched  off. 

The  previous  formulation,  i.e.  current  -i-  magnetic  field  on,  will  be  denoted  as  Case  0.  In  all  three 
cases  the  computations  were  performed  for  45mks  (VOmks  after  initiation).  The  two-dimensional 
fields  of  temperature  and  water  concentration  are  presented  in  Fig.4.  Initial  state  corresponding  to 
r=25mks  is  shown  in  Fig.4a.  In  Fig.4b  the  fields  at  r=70mks  are  presented  for  Case  0. 

Fig.4c  corresponds  to  the  Case  1  and  the  same  time  t=70mks.  Fig. Id  shows  the  fields 
corresponding  to  Case  2  at  t=70mks.  In  all  pictures  the  temperature  ranges  from  rmin=300K  to 
Tmax=4300K,  A7’=50K.  The  water  mole  fraction  varies  from  Cmm=0  to  Cmax=0.24,  AC=0.005. 

In  the  Case  0  the  largest  combustion  zone  is  observed.  The  arc  motion  increases  mixing 
which  provides  more  intensive  combustion.  In  the  Case  1  the  arc  slows  down  and  mixing  is  mainly 
due  to  diffusion  of  species.  In  the  Case  2  the  arc  kinematics  is  nearly  same  as  in  Case  1.  The  size  of 
the  arc  is  slightly  larger  in  Case  1  than  in  Case  0  because  of  Joule  heating.  Maximum  temperature  is 
nearly  same  as  in  Case  0,  4000K.  In  Case  2  the  maximum  temperature  is  about  1700K.  However 
this  large  difference  doesn’t  play  any  role,  since  combustion  initiates  at  rather  moderate 
temperatures.  One  can  say  that  in  no  current  case  the  temperature  is  maintained  only  due  to 
combustion.  The  total  production  of  water,  which  is  considered  as  the  combustion  efficiency,  is 
approximately  equal  for  both  no-interaction  cases.  This  is  well  seen  from  Fig.5,  where  the  total 
water  production  is  given  against  the  time.  The  Curve  1  on  the  plot  corresponds  to  Case  0,  Curve  2 
corresponds  to  Case  1,  and  Curve  3  corresponds  to  Case  2.  Maximum  production  is  observed  for 
Case  0  because  of  more  developed  interface  surface  due  to  arc  motion.  Another  feature  of 
combustion  process  can  be  seen  from  Fig.3,  where  the  total  enthalpy  of  the  fluid  volume  vs.  time  is 
shown  for  the  same  cases  as 
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Fig.4.  Temperature(left)  and  water  concentration(right)  for  initial 
state(a);  Case  0(b),  Casel,  Case  2;  t=70mks 
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those  in  Fig.5.  While  energy  input  in  Case  1  increases  in  time,  the  water  production,  as  combustion 
efficiency  sensor,  is  not  effected  by  this  circumstance.  The  maximum  production  rates  are  observed 
at  rather  moderate  temperatures.  At  high  temperatures,  say  above  2500K,  dissociation  processes 
prevail. 


Fig.5.  Total  water  production  for  Case  0(Curve  1),  Case  1  (Curve  2)  and  Case  2(Curve  3). 

Thus,  one  could  conclude  that  MHD-interaction  does  effect  on  mixing  and  combustion.  At 
first  turn,  this  is  due  to  the  motion  of  the  arc.  Electromagnetic  forces  increase  the  total  vorticity  in 
fluid  and  make  mixing  more  intensive.  The  heat  input  seems  to  play  the  role  of  starter:  mixing  of 
components  should  be  performed  at  even  moderate  temperatures  (above  1200K);  then  self-sustained 
combustion  could  start.  Initial  mixing  and  igniting  of  mixture  can  be  effectively  performed  with 
MHD-interactions . 
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Fig.6.  Total  enthalpy  for  Case  0  (Curve  1),  Case  1  (Curve  2)  and  Case  2  (Curve  3). 
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IV-3.  Numerical  prediction  of  arc-driven  mixing  and  combustion  under  experimental 
conditions 

The  previous  work  has  been  intended  to  study  the  capabilities  of  the  arc  discharge  driven  by 
the  magnetic  field  to  enhance  the  mixing  and  combustion  in  the  non-premixed  mixtures.  The 
simplified  problem  has  been  considered,  in  which  hydrogen  (fuel)  and  air  (oxidizer)  were  originally 
separated.  Mixing  and  combustion  were  stimulated  by  the  electric  arc  moved  by  the  applied 
magnetic  field  across  the  fuel-oxidizer  interface.  It  was  numerically  shown  that  the  motion  of  the 
arc  increases  the  contact  surface  while  the  high  temperature  within  arc  provides  ignition  and 
combustion  around  the  interface.  It  was  found  that  combustion  occurs  during  rather  long  time  after 
switching  off  the  electric  current. 

These  results  stimulated  us  to  perform  the  next  series  of  numerical  experiments.  The  goal  is 
to  study  combustion  in  non-premixed  fuel-oxidizer  jets.  The  numerical  formulation  is  assumed  to  be 
close  the  conditions  of  small-facility  experiments.  Such  facility  is  now  being  designed.  The 
interaction  of  two  flat  jets  is  considered.  The  electric  arc  normal  to  the  flow  is  generated  by  the 
electrode  system  so  that  the  arc  current  is  constant.  Magnetic  field  created  by  the  coil  is  directed 
along  the  flow.  Such  configuration  enforces  the  arc  to  move  across  the  jet  interface.  In  numerical 
predictions  we  shall  consider  the  following  formulation.  Let  two  jets  flow  in  X-direction.  The  jet  of 
hydrogen  is  initially  located  in  the  lower  half-plane,  and  the  jet  of  air  occupies  the  upper  half-plane. 
The  size  of  the  domain  in  X-direction  is  taken  to  be  4cm,  the  size  of  each  jet  in  Y-direction  is  taken 
to  be  0.5  cm.  The  flat  velocity  profiles  are  assumed  for  simplicity  at  the  inlet  station.  On  the  rest  of 
the  domain  boundary  the  outflow  conditions  are  specified.  Initial  velocity  field  corresponds  to  the 
inlet  conditions  overall  the  domain.  Z-directed  electric  arc  of  2mm  in  diameter  is  initiated  at  3  mm 
downstream  the  inlet  station  and  1  mm  below  the  jet  interface.  Initial  temperature  within  the  arc 
was  taken  4000K,  which  is  not  so  important  as  found  from  previous  calculations.  Duration  of  the 
arc  discharge  was  taken  20  microseconds,  and  constant  total  current  lOA  was  set.  The  constant  X- 
directed  magnetic  field  of  1  Teslas  was  assumed.  The  following  flow  parameters  were  used: 
pressure  -  Ibar,  temperature  -  300K  for  both  jets.  Two  sets  of  jet  velocities  were  considered. 

Case  1:  both  air  and  hydrogen  inlet  velocities  were  set  to  100  m/s; 

Case  2:  air  inlet  velocity  equals  to  100  m/s  and  hydrogen  inlet  velocity  equals  to  200m/s. 
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Computational  Details 


The  problem  formulated  above  was  solved  by  the  control-volume  method.  Exact  Riemann 
solver  was  used  to  evaluate  convective  fluxes  at  the  cell  boundaries.  Diffusive  fluxes  were 
estimated  by  the  method  similar  to  central-difference  scheme  on  the  orthogonal  grids.  Four-stage 
Runge-Kutta  time-stepping  was  applied  with  constant  CFL-number  equal  to  2.  The  chemical 
source-terms  (rates  of  species  production)  were  evaluated  in  implicit  manner  using  special  kinetic 
solver.  The  computational  grid  consisted  of  200  cells  in  flow  direction  and  100  cells  in  transversal 
one. 
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Figs  7.  Temperature  and  velocity  vectors,  Case  1. 
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Figs  8.  Water  mass  fraction  distributions.  Case  1. 

The  main  features  of  flow  can  be  seen  from  Fig.  7  and  Fig.  8  corresponded  to  Cases  1  and  2, 
respectively.  Fig.7  shows  temperature  and  velocity  fields  for  Case  1  at  different  time  moments. 
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Fig.8  shows  the  distribution  of  water  mass-fraction  at  the  same  time  moments.  Fig.9  and  Fig.  10 
represent  the  analogous  distributions  for  Case  2.  The  time  values  are  given  in  the  pictures. 
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Figs.  9.  Temperature  and  velocity  vectors,  Case  2. 
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Fig.lO.  Water  mass  fraction  distributions,  Case  2. 

The  characteristic  feature  of  the  flows  under  consideration  is  that  originally  small  arc  creates 
the  spot  with  well-developed  reaction  zone,  in  which  combustion  can  support  the  existence  of  the 
spot  during  a  long  time  compared  to  the  arc  life-time. 
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The  formulations  considered  above  were  used  as  the  first-stage  predictions  to  parameters  of 
the  small-scale  facility  being  designed.  In  a  future  coupled  experimental  work  and  numerical 
analysis  is  hoped  to  give  more  knowledge  about  the  topics  of  interest. 
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IV-4.  The  Reacting  Volume  and  Vorticity  Evolution  due  to  the  Body-Force 

In  order  to  prove  the  reacting  volume  concept  used  in  this  work  for  all  problem  formulations 
let  us  consider  the  quantitative  description  of  the  process  in  question. 

The  reacting  volume  value  can  be  define  as 

K=II,Xi.o'LXj,FdV; 

Vo  i  J 

where  and  Xj^p  are  the  mole  fraction  of  i-th  oxygen-containing  species  and  the  mole  fraction  of 
j-th  fuel-containing  species.  The  summing  should  be  made  on  all  species  in  the  mixture.  In 
principle,  such  a  definition  of  the  reacting  volume  value  could  be  normalized  with  a  limiting  value 
corresponding  to  uniform  mixture.  However,  the  constant  factor  omitted  in  the  expression  written 
above  is  non-important  in  the  further  analysis.  One  can  see  also  that  at  the  initial  state  when  two 
uniform  volumes  of  fuel  and  oxidizer  are  put  into  contact  the  reacting  volume  Vr  defined  here  is 
equal  to  zero.  It  is  also  clear  that  the  reacting  volume  is  a  measure  of  “combustion”  reactions  rate 
between  fuel-  and  oxygen-containing  molecules. 

Two  others  parameters  will  be  used  further  in  this  section. 

The  first  one  is  a  species  production 

e,=  JC,JV; 

where  C,  is  the  weight  fraction  of  i-th  species.  In  this  particular  case  the  two  product  species  will  be 
used:  1^0  and  OH,  to  characterize  the  overall  combustion  productivity.  Again,  non-normalized 
values  are  acceptable  for  the  particular  purposes  of  the  analysis. 

The  characteristics  of  the  flow  field  disturbances  created  by  a  body-force  jxB  is  given  by 
the  integral  vorticity  defined  as 

0=  J|rotv|^^y. 

This  parameter  represents  the  integral  vorticity  of  the  flow  field  in  the  domain  that  is 


responsible  for  the  enhancement  of  the  kinematics  mixing  due  to  the  stretching  of  a  fluid  volume  of 
the  media.  The  mechanism  of  this  process  is  described  in  details  in  our  preceding  work  [1]. 
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As  an  example  of  MHD-assisted  mixing  and  combustion  the  computational  domain  and  the 
initial  and  boundary  conditions  similar  to  those  considered  in  the  Sub-Sections  IV-1  and  IV-2  are 
utilized. 


Fig.  11.  Flow  Field  and  Fuel  Mole  Fraction  Distribution  at  t  =  70  mks 


In  the  Fig.  1 1  the  flow  field  structure  (the  right  half  of  the  domain  is  shown  because  of  the 
symmetry)  and  the  fuel  mole  fraction  (hydrogen)  are  presented  at  the  time  moment  t=70psec.  It  is 
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clear  from  the  picture  that  the  reacting  volume  represented  here  by  the  colored  area  between  the  red 
(100%  fuel)  and  blue  (0%  fuel)  is  increasing  significantly  due  to  the  vortex  flow  field  ereated  by  the 
body  force.  In  Fig.  12  the  eorresponding  value  of  the  reacting  volume  value  Vr  defined  above  is 
shown.  The  irregular  behavior  at  the  initial  stage  is  caused  probably  by  the  starting  shock  wave 
created  by  the  heat  release  in  the  arc  and  the  following  the  mass  loss  through  the  domain  border. 
After  that  the  practically  linear  growth  of  the  reacting  volume  is  observed.  It  is  interesting  that  total 
mass  production  of  the  eombustion  ignited  by  the  heating  is  also  praetically  linear  in  respect  to 
time.  One  can  conclude  that  the  mixing  process  (the  formation  of  the  reacting  volume)  is  of  primary 
importance  for  the  eombustion  efficieney.  At  the  same  time  it  should  be  noted  the  1^0  production 
reveals  a  slight  accelerating  due  to  probably  cooling  at  the  periphery  of  the  reacting  area  providing 
the  conversion  of  OH  to  H2O. 

The  information  supporting  this  conclusion  is  presented  in  Fig.  13  where  the  vorticity 
generation  is  plotted.  The  irregular  interval  at  the  beginning  of  mixing  and  combustion  induced  by 
MHD  body  force  is  even  more  elearly  seen  from  this  plot. 


Fig.  12.  Reacting  Volume  and  Mass  production  of  1^0  and  OH  vs  time 
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Fig.  13.  Vorticity  Evolution  in  Time 


Fig.l4.  Reacting  volume  vs.  time.  Curve  1 :  diffusion  only,  constant  temperature  300K;  Curve  2: 
diffusion  only,  hydrogen  temperature  3000K;  Curve  3:  diffusion  and  combustion,  hydrogen 
temperature  300K;  Curve  4:  MHD-assisted  combustion. 
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To  light  up  the  reacting  volume  concept  let  us  consider  the  simple  test.  The  same  initial 
fuel/oxidizer  state  is  considered,  i.e.  hydrogen  is  “specified”  in  the  bottom  part  of  the  domain  and 
the  oxidizer  (oxygen)  is  specified  in  the  top  part.  The  system  is  maintained  at  the  room  temperature. 
The  reacting  volume  in  this  case  is  the  potential  zone  where  the  combustion  could  take  place  if  it 
were  ignited.  The  reacting  volume  is  exactly  the  diffusion  layer  which  is  small  and  grows  in  time  as 
square  root.  The  development  of  no-combustion  reacting  volume  is  presented  in  Fig.  14  by  the 
Curve  1.  The  temperature  dependence  of  the  no-combustion  reacting  volume  is  demonstrated  by  the 
Curve  2  in  the  same  plot.  This  curve  stands  for  the  initial  temperature  of  hydrogen  3000K.  The 
larger  reacting  volume  as  compared  to  the  previous  case  is  due  to  enhanced  diffusion  coefficients 
because  of  higher  initial  temperature.  The  Curve  3  in  Fig.  14  corresponds  to  the  initial  state  of  the 
previous  case,  but  with  the  combustion  allowed.  Slight  increase  of  the  reacting  volume  id  thought  to 
be  due  to  expansion  because  of  heat  release.  For  comparison,  the  reacting  volume  corresponding  to 
the  MHD-assisted  combustion  from  the  previous  graph  is  also  shown  as  Curve  4.  The  difference 
between  MHD-assisted  combustion  and  non-MHD  one  is  evident. 
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V.  EXPERIMENTAL  FACILITY  PRELIMINARY  DESIGN 


V-1.  Experimental  Scheme  Development 

The  experimental  validation  of  the  proposed  MHD  assisted  mixing  mechanism  will  be  a 
subject  of  the  future  work  of  the  year  activity  of  the  Project.  At  this  stage  a  preliminary  analysis 
of  the  possible  approaches  is  undertaken. 

The  schematic  diagram  oft  he  proposed  experimental  setup  is  presented  in  Fig.l  1. 

The  main  components  are  as  following: 

1 .  High  pressure  vessels  for  working  gas  storage; 

2.  The  nozzles  for  formation  of  two  co-flow  streams  of  the  gas  to  be  mixed; 

3.  Electrodes  system  to  provide  pulse  electrical  discharge  in  the  flow; 

4.  Pulse  magnetic  system  to  create  magnetic  field  in  proper  direction  in  the  mixing  volume; 

5.  Power  supply  systems  of  electrical  discharge  and  the  magnetic  system; 

6.  Diagnostics/visualization  equipment; 

7.  Ventilation 

It  is  proposed  to  use  at  the  first  stage  two  contrast  colored  gases  such  an  air  and  smoked  gas. 
Another  possibility  acceptable  in  a  case  of  utilization  of  the  Schlieren  technique  could  be  a  mixing 
experiment  with  two  streams  with  different  temperature  (density). 

The  proposed  scheme  of  the  experiment  consists  of  the  several  steps: 

1.  The  initialization  of  the  co-flow  free  jets  in  the  test  volume; 

2.  The  initialization  of  the  pulse  discharge  across  the  streams  parallel  to  the  interface  between 
two  gases; 

3.  The  synchronized  induced  magnetic  field  should  be  oriented  perpendicular  to  the  main 
discharge  current,  and  correspondingly  parallel  to  the  interface; 

4.  The  visualization  system  should  record  the  effect  of  the  j*B  pondermotive  force  acting  in 
the  vicinity  of  the  interface. 

At  the  next  stage  the  ignition  and  combustion  process  affected  by  the  pondermotive  force 
and  concentrated  heat  release  due  to  Joule  dissipation  will  be  studied  as  well. 
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The  experimental  parameters  are  estimated  with  utilization  of  the  results  of  the  theoretical 
analysis  and  numerical  simulation. 

As  it  was  already  mentioned  in  the  previous  Interim  report  the  typical  value  of  the  current 
should  be  around  10  Amps,  and  the  magnetic  induction  is  about  2-3T.  These  parameters  provide  the 
current  channel  displacement  in  several  millimeters  resulting  in  penetration  of  heated  portion  of  one 
component  to  another  co-flow  stream  during  the  observation  time.  The  latter  is  defined  by  the 
velocity  of  the  flow  and  the  size  of  observed  area. 

The  discharge  current  will  be  provided  by  high  voltage  (30-50  kV)  power  supply  -  the 
standard  high  voltage  transformer  with  rectifier  BT-140.  This  transformer  will  be  used  to  charge 
high  voltage  capacitor  ?  ?  ?  -12  of  12  mkF  through  two  load  resistors  ???-20  of  4.3  MOhm  each. 
This  could  provide  a  discharge  current  pulse  about  100 Amps,  at  voltage  30  kV  and  duration  about  1 
msec.  The  estimated  interelectrode  gap  is  about  3  cm. 

The  magnetic  system  will  consists  of  the  coil  of  5  turns  with  6  cm  in  diameter.  The  coil 
current  of  10  kAmps  will  result  in  the  IT  magnetic  induction  in  the  test  volume.  The  power  supply 
of  such  a  magnet  will  be  solved  with  a  standard  high  voltage  power  supply  of  BT-20  that  will 
charge  a  capacitor  bank  of  20  K417-7  units  with  total  capacity  of  2  mF  at  the  voltage  5  kV.  In  order 
to  provide  the  pulse  duration  more  than  1  msec  (that  is  a  synchronization  condition  with  discharge 
current  duration)  the  total  inductance  should  be  not  less  than  12  mkH.  Thus,  the  inductance  of 
system  can  be  fitted  with  the  coil  turn  number  used.  It  should  be  noted  that  increasing  the  number  of 
coil  turn  results  in  increasing  of  magnetic  induction  that  is  very  desirable  for  the  final  purposes  of 
the  experiment. 

The  synchronization  of  the  test  consequence  will  be  provided  with  the  CAMAC  system 
available  at  the  site. 
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Fig.4.  Layout  of  experimental  rig  on  MHD  assisted  mixing 
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V-2.  Experimental  rig  assemblage. 


The  schematic  diagram  of  preliminary  proposed  experimental  facility  presented  and  briefly 
discussed  in  the  previous  Sub-Section  is  used  now  in  preparation  of  the  experimental  study  of  the 
MHD  driven  mixing  and,  the  second,  stage  ignition  and  combustion.  At  the  present  time  the 
experimental  facility  is  being  to  be  assemblage  and  tuned.  The  first  stage  is  aimed  to  create  the 
simplest  working  version  of  the  experiments  to  check  the  different  options  of  the  experimental 
program  and  arrangement. 

In  Photos  3.1  through  3.3  the  current  view  of  the  assemblage  is  presented. 


Fig.3.1  The  Experimental  Assemblage  of  the  MHD  Assisted  Mixing  Study 


The  main  components  are  as  follows.  1  -  High  voltage  supply  BC20  for  the  voltage  0-20  kV;  2  - 
capacitor  of  lOOmkF  (two  piecesO;  3  -  high  voltage  voltmeter;  4  -  quartz  tube  of  exhausting 
system;  5  -  example  of  the  magnetic  coil  strength  elements;  6  -  the  fast  camera  ????-16  of 
5000fr-sec  maximal  n  7  -  magnetic  system  installed;  8  -  mercury  switcher  ???-6;  9  -  converging 
nozzle;  10  -  ventilator;  1 1  -  the  fast  camera  controller. 
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Fig.3.2.  The  side  view  of  the  experimental  set-up 
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Fig.3.3  The  nozzle  and  the  magnetic  coil  view  in  assemblage 


The  components  of  the  experimental  set-up  are  tested  now. 
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CONCLUSION 


The  new  concept  of  MHD  driven  Advanced  mixing  of  gaseous  fuel-oxidant  system  in  a 
high-speed  combustor  is  formulated. 

The  background  of  such  a  concept  consists  of  the  earlier  developed  model  of  kinematics 
mixing  in  vortex  turbulent  flows. 

The  vortex  nature  of  the  Lorentz  body  force  resulting  from  perpendicular  vectors  of 
electrical  current  and  externally  applied  magnetic  filed  provides  the  flow  vorticity  intensification. 

The  enhancement  of  the  mixing  under  such  conditions  is  characterized  by  an  exponential 
growth  of  the  contact  surface  square  in  a  turbulent  flow. 

The  second  mechanism  -  molecular  diffusion  is  responsible  for  the  creating  the  flammable 
mixture  in  a  physically  small  volume. 

These  two  mixing  mechanisms  are  combined  to  formulate  the  concept  of  reacting  volume  - 
the  volume  where  the  combustion  is  only  possible. 

Besides  of  the  qualitative  physical  model  the  proper  numerical  model  has  been  developed  to 
study  the  process  on  the  microscopic  scale. 

This  model  includes  the  treatment  of  chemical  reacting  flows  on  the  base  of  full  Navier- 
Stockes  approach,  chemical  kinetics  model,  and  MHD  effects  model. 

Preliminary  results  confirm  the  basic  physical  ideas.  The  combustion  in  co-flow  streams  can 
be  effectively  initiated  and  controlled  by  the  manipulation  of  the  electrical  parameters:  strength  of 
external  magnetic  field  and  electrical  discharge  current. 

The  dramatic  growth  of  the  reacting  volume  and,  hence,  the  combustion  performance  in  the 
presence  of  MHD  interaction  was  clearly  demonstrated. 

The  preliminary  sketch  of  the  future  laboratory  scale  experiment  is  presented  as  well. 

The  parametric  study  of  the  proposed  concept  as  applied  to  the  conditions  typical  for  the 
practical  situation  in  a  scramjet/  ramjet  will  be  completed  in  a  related  project  supported  by  ISTC  in 
2000/2001. 
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Nomenclature 


J 

B 

L,l 

t 

U 


F,G 

Fv,G, 

S 


P 

y=iyx,Vy) 

E 

Yi 

Xi 

P 


3i~(JxijJyi) 


X 

T 

hi 

kf 

Dim 

N 


Da 


Mi 


W,M 


-  electric  current  density 

-  magnetic  induction  vector 

-  characteristic  lengths 

-  time,  characteristic  time 

-  conservative-variable  vector  consisting  of  fluid  density,  momentum,  total  energy  and 
species  densities 

-  inviscid  flux- vectors  in  x-  and  y-directions,  respectively 

-  viscous  flux- vector  in  x-  and  y-directions 

-  source-term  vector 

-  fluid  density 

-  fluid  velocity 

-  total  specific  energy 

-  species  mass-fractions 

-  species  mole-fractions 

-  thermodynamic  pressure 

-  viscous  stress-tensor  components 

-  mass  diffusion  flux  of  species  i 

-  energy  flux  due  to  heat  conductance,  viscous  dissipation  and  multi-component 
diffusion 

-  molecular  viscosity  coefficient 

-  molecular  heat  conductivity  coefficient 

-  temperature 

-  specific  enthalpy  of  species  i 

-  formation  enthalpy  of  species  i 

-  effective  mass-diffusion  coefficient  of  species  i 

-  number  of  species 

-  binary  diffusion  coefficients 

-  molar  weight  of  species  i 

-  molar  weight  of  mixture 


o. 


-  characteristic  cross-section  of  i-j  pair  collision, 


J 


^G-i)  _  i_j  collision-integral  of  (1.1)  order  (viscous  type) 

T*  -characteristic  temperature  in  the  Lennard-Jones  inter-molecular  interaction  potential 

^C2.2)  _  collision-integral  of  diffusion-type 

pP  -  universal  gas  constant 

Ri=lP/Wi  -  gas  constant  of  species  i 

Cpi  -  specific  heat  at  constant  pressure  of  /-th  species 

Cvi  -  specific  heat  at  constant  volume  of  /-th  species 

N 

R=^  YjRj  -  gas  constant  of  mixture 
i=l 

N 

h  =  -  specific  enthalpy  of  mixture 

i=l 

e  =  h-  pjp  -  specific  internal  energy 

P  -  Hall  parameter,  P=a)eX^,  (O  ^  is  the  gyro-frequency  =  eBlnie,  x  ^  is  the  average  time 

between  impacts  for  electron 
a  -  electrical  conductivity 

E  -  electric  field  strength  vector 

Sx,y  -  components  of  electromagnetic  force  [jxB] 

Se  -  energy  source  due  to  work  of  electromagnetic  force  and  radiation 

Qrad  -  radiation  power 

<3 SB  -  Stefan- Boltzman  constant 

y"i,r  -  the  stoichiometric  coefficients  of  reagents  and  products,  respectively 
Ci  -  molar  concentrations  of  species  /,  c,=p  y;/W, 

kfr ,  kbr  -  forward  and  backward  reaction  rate  constants 
Wj  -  /-th  species  rate  of  mass  due  to  chemical  reactions 

CO  -  computational  domain  volume 

9(0  -  boundary  of  computational  domain 

R  -  residual-vector  in  the  discrete  Navier-Stokes  equations 

r”  -  time  value  on  n-th  time-integration  step 

-  time-integration  step 

U",  -  solution-vector  at  t"  and  time-level,  respectively 

Ui  -  intermediate  solution-vectors  in  Runge-Kutta  method 

Ul  ,  Ui;  -  left  and  right  solution- vectors  respective  to  integration-point  (ip)  position 


\|/p,£  -  limiter-functions  in  linear-in-cell  variable  interpolation 

larc  -  total  arc  current,  discharge  current 

Vr  -  reacting  volume  measure 

Qi  -  total  production  of  species  i 

Q  -  integral  vorticity  of  the  flow 


